Charged Impurity Scattering in Graphene Nanostructures 
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We study charged impurity scattering and static screening in a top-gated substrate-supported 
graphene nanostructure. Our model describes how boundary conditions can be incorporated into 
scattering, sheds light on the dielectric response of these nanostructures, provides insights into the 
effect of the top gate on impurity scattering, and predicts that the carrier mobility in such graphene 
heterostructures decreases with increasing top dielectric thickness and higher carrier density. An 
increase of up to almost 60 percent in carrier mobility in ultrathin top-gated graphene is predicted. 



Introduction- The electron mobility in single-layer 
graphene (SLG) has been demonstrated to be as high 
as 200,000 cm^V~"'^s~^ and is central to many of its 
potential nanoelectronic applications. However, in most 
graphene-based heterostructures, SLG must be physi- 
cally supported by an insulating dielectric substrate such 
as Si02, and the carrier mobility in such structures is 
about one order of magnitude lower |lj] as a result of scat- 
tering with impurities, defects, and surface roughness. 
This reduction in carrier mobility can be further exac- 
erbated when a top gate, consisting of a layer of high-K 
dielectric material such as Hf02 or AI2O3 overlayed with 
metal, is deposited on SLG 0-131 ■ 

Experiments suggest that the dominant factor limiting 
electrical transport is scattering by charged impurities, 
which are believed to be located at or near the SLG- 
substrate interface [sl]. Adam and co-workers have suc- 
cessfully explained the linearity of the conductivity with 
respect to carrier density as a consequence of charged 
impurity scattering 0, 0|- However, their theory is lim- 
ited to the simple geometry of SLG supported by a sub- 
strate. In more realistic graphene heterostructures, the 
SLG is encapsulated between the substrate and the top 
gate d, H, Q , an arrangement which offers better local 
electrostatic control than the bottom gate and is prob- 
ably required for large scale integration. On the other 
hand, top gates modify the dielectric environment of the 
SLG by introducing image charges along the various in- 
terfaces. Thus, given the significance of impurity scat- 
tering, it is important to treat this process theoretically 
by accounting fully for these geometrical effects in order 
to optimize the design of SLG-based nanoelectronics. 

Jena and Konar suggest that charged impurity 
scattering in low-dimensional semiconductor nanostruc- 
tures can be damped by coatingthem with high-K di- 
electrics. It has been shown [111 . 1 131 that placing a 
higher- K overlayer {e.g., ice in Ref. |ll| . non-polar liq- 
uid in Ref. [12]) on SLG can lead to a weakening of the 
Coulombic scattering forces and can lead to an improve- 
ment in the conductivity (although Ponomarenko and co- 
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workers using a high-K liquid overlayer were not able find 
any significant improvement fis']). The deposition of di- 
electric films, such as Hf02 or AI2O3, on graphene has 
been found to lead to the degradation of electron mobil- 
ity [3, [l3| . This has been attributed to the roughening 
of the graphene surface. However, using atomic layer de- 
position (ALD), Fallahazad and co-workers Q have 
been able to deposit ultrathin high-K dielectric materials 
on SLG with much less roughness and observe a signifi- 
cant improvement in the carrier mobility. Hollander and 
co-workers were also able to see an increase of the Hall 
mobility in epitaxial graphene with thinner top gate di- 
electrics (l5| . 

In this paper, we extend the theory of Adam and co- 
workers to the more complicated geometry of top-gated 
SLG, and apply it to study the effect of the top gate on 
impurity scattering, emphasizing the dependence of the 
mobility on the top dielectric thickness. Our principal 
finding is that a thinner top gate leads to a higher car- 
rier mobility because of the stronger image-charge effect 
in the gate. This agrees with the scaling trend found 
by Fallahazad and co-workers [sll as well as Hollander 
and co-workers [151]. This mobility improvement due to 
stronger gate screening was also suggested for conven- 
tional metal-oxide-semiconductor systems by Gamiz and 
Fischetti [l6j . A general formula for charged impurity 
scattering in SLG heterostructures is provided. 

The outline of the paper is as follows. We derive 
the Thomas-Fermi screening and the electron-impurity 
Coulombic interaction, taking into account the bound- 
ary conditions imposed by the substrate and the top gate. 
The matrix element for electron-impurity interaction is 
calculated and used to compute the momentum relax- 
ation rate which determines the conductivity and carrier 
mobility. We then compute the carrier mobility for dif- 
ferent top dielectrics, gate thickness, and carrier density. 

Method- The model of the double-gated graphene 
nanostructure consists of a SLG sheet sandwiched be- 
tween two oxide layers with the interfaces coplanar with 
the x-y plane. The substrate consists of the semi-infinite 
region z < while the top oxide spans the region 
h < z < h + tox , where h is the height of the empty space 
between the two oxide layers and tox is the thickness of 
the top oxide layer. The region z > h + tox is assumed 
to be composed of a semi-infinite ideal metal. Thus, all 
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Figure 1: Basic model used in our calculation. The SLG is an 
infinitely thin layer suspended in the middle of the vacuum 
gap between a semi-infinite substrate and a top oxide layer of 
thickness tox- The dielectric is capped with metal which we 
assume to be a perfect conductor. The charged impurity on 
the SLG has image charges under and above it in the substrate 
and top gate respectively. 



field lines terminate at z = h + tox- For simplicity, we 
assume that the SLG floats at a height oi z = d = h/2 
halfway between the oxide layers. The structure is shown 
in Fig. [1] The dielectric constants of the top and bottom 
oxides are etox and £box respectively. 

To describe the static dielectric screening of a charged 
impurity in SLG, we start from the Poisson equation for 
the screened scalar potential ^scr'- 
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.(R,z)] , (1) 



where pimp is the impurity charge and pscr represents 
the screening (polarization) charge. Here R represents 
the coordinate on the x-y plane while z is the coordinate 
along the perpendicular direction. The integral form of 
Eq. HI is: 

$,cr(R,2) = <I>(R,z) + J dR'dz'G{Rz,R'z')p,cr{R'z') , 

(2) 

where G{Rz, R'z') is the Green function that satis- 
fies the equation -V^ [e(R, z)G'(R2, RV)] = S{R ~ 
R' , z — z') with the proper boundary conditions. 
The bare potential $(R, z) is defined as $(R, z) — 
J dR'dz'G{Rz,R' z')pimp{R' z'). The second term on the 
right-hand side of Eq. ([2]) represents the screening charge 
distribution. The bare and screened potentials can be 
written in terms of their Fourier components: $(R, z) = 
Y.ciMzy-''^-^ and ^scr{R.z) = EQ<^Q'^(^)e-*'^-^. 
The two-dimensional Fourier transform of Eq. ^ yields 
the following expression for the z-dependent part of the 
Fourier-transformed screened potential: 



'/'q''(^) = '^q(^) + j dz'GQ{z,z')p^^-{z') , (3) 

where Q is the two-dimensional in-plane wave vector. 
Equation ^ can be solved if we express the polarization 



charge pg'" as a function of the screened scalar potential. 
Thus, we express the screening charge as: 



p^^'-(z) = e2n(g)/(z)0^-(z) , 



(4) 



where Ii{Q) is the in-plane static polarizability and f{z) 
governs the polarization charge distribution in the per- 
pendicular direction. We require f{z) to be normalizable 
to unity, i.e., / dz f{z) — 1. For convenience, we choose 
f{z) — 6{z — d). Physically, this implies that the SLG is 
idealized as an infinitely thin sheet of polarized charge. 
Combining Eqs. ([3]) and we obtain the expression: 

't^Q^'iz) = 0q(^) + 6^ / dz'GQ{z,z')n{Q)f{z'W^^{z') . 

(5) 



After some algebra, Eq. ^ yields: 

e^GQiz,d)Il{Q) 



s2GQ(d,d)n(Q) 



Since we are only concerned about the in-plane scattering 
potential (j)Q'^{d), we have: 
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l-e^Ggid, d)U{Q) 



(6) 



Therefore, the generalized expression for the static 
two-dimensional screening function is e{Q) = 1 — 
e^GQ{d, d)Il{Q). In the simplest case of suspended SLG, 
we have Ggid^d) = ^ and n(g) = giv- 
ing us e(Q) = 1 + ^, where = T^M^^ 
vp are the inverse screening length, the Fermi level and 
the Fermi velocity, respectively. Note that the screening 
function is determined by its electrostatic environment 
through GQ{d,d), the electrostatic Green function. 

In order to solve Eq. ([5]) and find the screening function 
e((5), we need an explicit expression for the electrostatic 
Green function Gq{z,z'). This Green function incorpo- 
rates the details of the electrostatic environment around 
the SLG. For a source in the empty space between the top 
and bottom oxide layers, i.e., for < z' < h, Gq{z,z') 
satisfies the equation: 

-e(^-g') Gq{z,z')^S{z~z') . 

This implies that we can always express Gq{z,z') as a 
linear combination of e"^^ and e~'^^ , or cosh((3z) and 
sinh(Qz). The Green function has to satisfy the follow- 
ing continuity conditions at the interfaces z = 0, z = /i 
and z = h + tox-Ggih + tox,z') = lim^^.oo Gq{z, z') = 
0, GQ{h+,z') - GQ{h-,z'), Gq(0+,z') = Gq(0-,z'), 

dGQ(z=h+,z') _ dGQ(z=h-,z') , aGQ{z=0+,z') _ 

etox dr~~ - ^0 55 '^0 Wz - 

gboa ^'^^^ In addition, we have the following con- 



ditions at z 
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Gq(2 



and 



eo^GQ(z = z'-,z') -eo^GQ(z - z'+,z') = 1. After 
some straightforward but lengthly algebra we obtain for 
Gq(0 < z < ft,, < z' < /i) the following expression: 



3 



Gq(z,z') = 



2eoQ 



1 - XtXbe~^Q'' 



'^-Q(z-z') _^^Qiz-z') 



X^eQ(z+z'-2h) ^ ^^^-Qiz+z') 

1 - Af A6e-2Q^ 



(7) 



where Xt = S ^"^j^Q*"--^" and = S^Z"" ■ In the 
region h < z < h + tox > making use of the continuity of 
the Green function and the fact that it has to terminate 
at z = /i + tox, we can easily write the Green function 
as Gq{z > h,z') = G{h,z' r''\'^^l^+l2''l Similarly, for 
z < we have Gq{z < 0, z') = G(0, zOe+Q^ 

The scattering potential due to a single bare point 
charge on the SLG at the origin satisfies the equation 
— eoV2$(R, d) = e^S{Il,d). The corresponding two- 
dimensional Fourier transform of $(R, d) is 4iQ{d) — 
e^GQ{d,d) . Upon including the screening function from 
Eq. we obtain the generalized expression for the 

screened Coulomb scatterer: 
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}{d,d) 



l-e^GQ{d, d)U{Q) 



(8) 



The matrix element of the scattering potential of ran- 
domly distributed screened impurity charge-centers in 



SLG is |(T4k,sK' 



l<^|K-K'|(^)l (1 + cosW)/2, 



where 0kk' is the angle between the wave vectors K and 
K', and s = +1 (-1) for the conduction (valence) band. 
The transport scattering rate is given by 1/t{Es-k) ~ 

{riimp/h) J dK'|(14K,.K')|'[l - COSflKK'I'JI-B.K - EsK'), 



where 



is the impurity density, and the conductivity 



a is approximated by cr = e'^VpD{Ep)T{Ep) /2, where 
D{Ep) is the density of states at the Fermi level. The 
impurity-limited mobility is taken to be 0, 0|: 
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(9) 



where n — Ep/{TTh'^Vp) is the zero-Kelvin approximation 
to the carrier density. 

Results and discussion- Using Eq. we estimate the 
impurity-limited mobility of top-gated SLG for different 
carrier densities, top gate dielectrics, and dielectric thick- 
ness. In particular, the dependence of the mobility on tox 
is studied in detail. For simplicity, we assume that the 
doping level is > 10^^ cm^^ because charge inhomogene- 
ity ('puddles') is a phenomenon at /near the Dirac point 
and is not captured in our theory. 

We consider three top gate dielectrics, Si02 (k=3.9), 
AI2O3 (k=12.5) dl and HfO2(K=22.0) @|. The bottom 
dielectric is assumed to be Si02. tox is varied from 1 to 
12 nm. The electron mobility is calculated for two values 
of the carrier density, 1 x 10^^ cm"^ and 5 x 10^^ cm^^. 
The results are shown in Fig. [2j We set the impurity den- 
sity to be Uimp = 0.5 X lO"'^^ cm"^. In Fig. [21 we notice 
that the dependence of the mobility on tox is stronger 
for n = 1 X 10^2 cm^2 than for n = 5 x 10^^ cm"^. As 
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Figure 2: Dependence of the carrier mobility (/i) on top gate 
dielectric thickness (tox) for Si02, AI2O3, and Hf02. The mo- 
bility decreases with increasing tox- This decrease is greater 
at n = 10^^ cm~^ than at n = 5 x 10^^ cm"^. 



expected, the higher-K oxides yield a larger mobility be- 
cause the higher effective permittivity weakens Coulom- 
bic interactions (electron-electron and electron-impurity) 
in graphene. Although the weaker electron-electron in- 
teraction results in less screening charge around the im- 
purity, the bare charge on the impurity is reduced to a 
greater extent, resulting in a overall diminished screened 
impurity charge. 

We also find that the mobility increases with decreas- 
ing tox, in qualitative agreement with the experimental 
results in Ref. [8]. In our simulations, this increase in 
mobility can be substantial. For example, in Fig. [2] at 
n = 10^2 cm^2, the mobility for a 1.0 nm Si02 top gate 
is almost 60 percent higher than that for a 12.0 nm top 
gate. This scaling trend has been attributed to the in- 
crease in defect concentration with greater tox [1| ■ How- 
ever, our model also reproduces this trend because with 
a decreasing tox the image 'anti-charge' at the metal- 
dielectric interface gets closer to the impurity charge 
and weakens its scattering potential. This implies that 
ultra-thin top gates can be used to compensate for large 
charged impurity densities. The carrier-density depen- 
dence of the mobility is calculated for different top gate 
dielectrics at tox= 2 nm and 12 nm. The carrier den- 
sity is varied from 10^^ to 10^'^ cm^'^. The results are 
shown in Fig. [31 In general, the mobility decreases with 
increasing carrier density, in contrast to the results in 
Refs. 01 which suggest that mobility should be inde- 
pendent of carrier density. This decrease is more acute 
in the 2 nm case than in the 12 nm case and for higher- 
K oxides. To understand why, we look at Eq. In 
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Figure 3: Dependence of the mobility on carrier density. The 
mobiUty decreases with carrier density, and the extent of the 
decrease is greater with smaller top gate dielectric thickness. 

suspended SLG, the Fermi level Ep is proportional to 
n^/^. The scattering time t(Ep) is also proportional to 
n^/^ because Gqid^d) in Eq. ([8]) scales as 1/Q. Thus, 
/i is independent of n. On the other hand, in top-gated 
SLG, Gq^d, d) is almost constant over a large range of Q 
values, and hence, t{Ef) is approximately independent 
of n at low carrier densities. Thus, /i scales roughly as 
71^^/^. Thus, the trend shown in Fig. [3]should be viewed 
as low-density/thin-oxide enhancement of the mobility 
above its gate-unscreened value thanks to the screening 
effect of the image charges induced in the gate by the 
impurities. 

The improvement of mobility with a thinner top oxide 



layer highlights the compensating effect of having a metal 
layer in close proximity to SLG. It has been suggested 
that voltage fluctuations induced by charged impurities 
in the substrate lead to local electrostatic doping and 
are responsible for the formation of 'puddles' in graphene 
near the charge neutral point (TtI . ITsI . We propose that 
doping inhomogeneities can be 'smoothened out' through 
the use of an ultrathin high-K top-gate structure since 
the image charges can partially neutralize the effects of 
the charged impurities. Unintentional doping from the 
substrate can also be reduced, leading to an improvement 
of carrier transport in top-gated SLG. 

Conclusions- We have studied charged impurity scat- 
tering and static screening in top-gated SLG nanostruc- 
tures. The image charge effect dampens Coulombic in- 
teraction, resulting in weaker unscreened charged impu- 
rity scattering and static screening. The carrier mobil- 
ity is found to increase with decreasing top gate dielec- 
tric thickness and decreasing carrier density. Our theory 
suggests that if charged impurity scattering is a signifi- 
cant factor in limiting electronic transport in SLG, then 
having a top metal gate with an ultrathin dielectric can 
compensate for the effect of the charged impurities. This 
offers a simple and practical design strategy to improve 
the impurity-limited mobility in SLG heterostructures. 
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